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Abstract: 

The cross section for the J/ ijj and T interaction with light hadrons is calculated in 
short-distance QCD, based on the large heavy quark mass and the resulting large energy 
gap to open charm or beauty. The low energy form of the cross section is determined 
by the gluon structure functions at large x; hence it remains very small until quite high 
energies. This behaviour is experimentally confirmed by charm photoproduction data. 
It is shown to exclude J/ ip absorption in confined hadronic matter of the size or density 
attainable in nuclear collisions; in contrast, the harder gluon spectrum in deconfined 
matter allows break-up interactions. 
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The interaction of heavy quark-antiquark resonances (quarkonia) with the ordinary 
light hadrons plays an important role both in the dynamics and in the thermodyna- 
mics of strong interaction physics. For QCD dynamics, it is important since the small 
quarkonium size probes the short-distance aspects of the big light hadrons and thus 
makes a parton-based calculation of the overall cross section possible [1-4]. In QCD 
thermodynamics, quarkonia can be used as a probe for deconfinement [5], provided their 
interaction in dense confined matter can be distinguished from that in a quark-gluon 
plasma. 

In this letter we will first recall the QCD analysis of quarkonium interactions with 
light hadrons. It concludes that the small size of quarkonia, combined with the rather 
large mass gap to open charm or beauty, strongly inhibits their break-up by low energy 
collisions with light hadrons [2]; the total quarkonium-hadron cross section attains a 
constant asymptotic value only at very high energies, compared to corresponding cross 
sections for the interaction between light hadrons. This slowly rising form of the cross 
section is derived from an operator product expansion with ensuing sum rules and 
becomes quite transparent in parton language. It is experimentally well supported by 
applying the corresponding analysis to data for charm photoproduction [6-9]. 

We then use this result to show that a slow J/ip or T cannot be absorbed (i.e., 
broken up into a DD) in hadronic matter of the volume and density which can be 
obtained in nuclear collisions. The essential point is here that absorption requires in- 
teractions with sufficiently hard gluons, and the gluon distribution within hadrons does 
not provide these for matter at meaningful temperatures. Hence absorption in confined 
hadronic matter is excluded as a possible cause of the J/ip suppression observed in 
nucleus-nucleus collisions. On the other hand, we find that quarkonia can be broken 
in deconfined matter, which contains much harder gluons. This provides a dynamical 
basis of quarkonium suppression by colour screening [5,10]. 

The QCD analysis of quarkonium interactions applies to heavy and strongly bound 
quark-antiquark states [2]; therefore we here restrict ourselves to the lowest cc and bb 
vector states J ftp and T, which we denote generically by following the notation of 
[2]. For such states, both the masses tuq of the constituent quarks and the binding 
energies eo($) — (2M(g g ) — M$) are much larger than the typical scale Aqcd for non- 
perturbative interactions; here (Qq) denotes the lowest open charm or beauty state. In 
$ — h interactions, as well as in $-photoproduction, 7/1 — > we thus only probe a 
small spatial region of the light hadron h; these processes are much like deep-inelastic 
lepton-hadron scattering, with large tuq and eo in place of the large virtual photon 
mass a/ —q 2 . As a result, the calculation of ^-photoproduction and of absorptive $ — h 
interactions can be carried out in the short-distance formalism of QCD. Just like deep- 
inelastic leptoproduction, these reactions probe the parton structure of the light hadron, 
and so the corresponding cross sections can be calculated in terms of parton interactions 
and structure functions. 

In the following, we shall first sketch the theoretical basis which allows quarkonium 
interactions with light hadrons to be treated by the same techniques as used in deep- 
inelastic lepton-hadron scattering or in the photoproduction of charm. We show the 
derivation for the sum rules which relate the absorptive $ — h cross section to hadronic 
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gluon structure functions [1,2]. This relation given, we calculate explicitly the energy 
dependence of the cross section. Readers only interested in this behaviour can therefore 
go immediately to eq. (24) . 

Consider the amplitude for forward scattering of a virtual photon on a nucleon, 

F(s, q 2 )~i J d 4 xe tqx < N\T{J^x)J„(Q)}\N > . (1) 

In the now standard application of QCD to deep-inelastic scattering one exploits the 
fact that at large spacelike photon momenta q the amplitude is dominated by small 
distances of order 1/a/~<7 2 (Fig- l a )- The Wilson operator product expansion then 
allows the evaluation of the amplitude at the unphysical point pq — » 0, where p is the 
four-momentum of the nucleon. Since the imaginary part of the amplitude (1) is propor- 
tional to the experimentally observed structure functions of deep-inelastic scattering, 
the use of dispersion relations relates the value of the amplitude at pq — > point to the 
integrals over the structure functions, leading to a set of dispersion sum rules [11]. The 
parton model can be considered then as a particularly useful approach satisfying these 
sum rules. 

In the case J M = Qj^Q, i.e., when vector electromagnetic current in eq. (1) is that 
of a heavy quark-antiquark pair, large momenta q are not needed to justify the use 
of perturbative methods. Even if q ~ 0, the small space-time scale of x is set by the 
mass of the charmed quark, and the characteristic distances which are important in the 
correlator (1) are of the order of l/2mg (Fig. lb). In [3,4], this observation was used 
to derive sum rules for charm photoproduction in a manner quite similar to that used 
for deep-inelastic scattering. 

In the interaction of quarkonium with light hadrons, again the small space scale 
is set by the mass of the heavy quark, and the characteristic distances involved are 
of the order of quarkonium size, i.e., smaller than the non-perturbative hadronic scale 
^qcd- Moreover, since heavy quarkonium and light hadrons do not have quarks in 
common, the only allowed exchanges are purely gluonic. However, the smallness of 
spatial size is not enough to justify the use of perturbative expansion [2]. Unlike in 
the case of ^-photoproduction, heavy quark lines now appear in the initial and final 
states (see Fig. lc), so that the (QQ) state can emit and absorb gluons at points 
along its world line widely separated in time. These gluons must be hard enough to 
interact with a compact colour singlet state (colour screening leads to a decoupling of 
soft gluons with the wavelengths larger than the size of the $); however, the interactions 
among the gluons can be soft and nonperturbative. We thus have to assure that the 
process is compact also in time. Since the absorption or emission of a gluon turns a 
colour singlet quarkonium state into a colour octet, the scale which regularizes the time 
correlation of such processes is by the quantum-mechanical uncertainty principle just 
the mass difference between the colour-octet and colour- singlet states of quarkonium: 
r c ~ l/(e8 — ei). The perturbative Coulomb-like piece of the heavy quark-antiquark 
interaction 

is attractive in the colour singlet (k = 1) and repulsive in the colour-octet (k = 8) state; 



in SU(N) gauge theory 

'/ JV 2 - 1 

« = <»> 
To leading order in 1/N, the mass gap between the singlet and octet states is therefore 
just the binding energy of the heavy quarkonium eo, and the characteristic correlation 
time for gluon absorption and emission is 

r c ~ l/e . (4) 

Although the charm quark is not heavy enough to ensure a pure Coulomb regime even 
for the lowest cc bound states (rj c and J/ip), the mass gap determined from the observed 
value of open charm threshold clearly shows that r c < Aqcd- For the T, the interaction 
is in fact essentially Coulomb-like and the mass gap to open beauty is even larger than 
for charm. One therefore expects to be able to treat quarkonium interactions with light 
hadrons by the same QCD methods that are used in deep-inelastic scattering and charm 
photoproduction. 

We thus use the operator product expansion to compute the amplitude of heavy 
quarkonium interaction with light hadrons, 

F$h = i I d 4 xe iqx (h\T{J(x)J(0)}\h) = ]T c n (Q, m Q ) (O n >, (5) 

J n 

where the set {O n } includes all local gauge invariant operators expressible in terms of 
gluon fields; the matrix elements (O n ) are taken between the initial and final light- 
hadron states. The coefficients c n are computable perturbatively [1] and process- 
independent. As noted above, in deep-inelastic scattering the expansion (5) is useful 
only in the vicinity of the point pq — > 0. The same is true for the case of quarkonium 
interaction with light hadrons. As shown in [2], the expansion (5) can therefore be 
rewritten as an expansion in the variable 

X= PQ_ = {s-Ml-Mp ^ (8 -Mp 

Mj> 2M$ 2M$ 1 ' 

where is the mass of the light hadron; the approximate equality becomes valid the 
heavy quark limit. For the lowest IS quarkonium state one then obtains 

oo / A \ n 

F* h = 4e 2 J2d n (O n )(-) , (7) 

where ro and eo are Bohr radius and binding energy of the quarkonium, and the sum 
runs over even values of n to ensure the crossing symmetry of the amplitude. The most 
important coefficients d n were computed in [1] to leading order in g 2 and 1/N. 

Since the total $ — h cross section a$>h is proportional to the imaginary part of the 
amplitude F$h, the dispersion integral over A leads to the sum rules 

- / dX \- n a* h (\) = r 3 e 2 d n (O n ) - . (8) 



Eq. (7) provides only the inelastic intermediate states in the unitarity relation, since 
direct elastic scattering leads to contributions of order 7q. Hence the total cross section 
in eq. (8) is due to absorptive interactions only [2], and the integration in eq. (8) starts 
at a lower limit Ao > Mh- Recalling now the expressions for radius and binding energy 
of IS Coulomb bound states of a heavy quark-antiquark pair, 

V 3# 2 / m Q 

eo = (m) mQ ' (10) 

and using the coefficients d n from [1] , it is possible [2] to rewrite these sum rules in the 
form 
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, (On), (11) 
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with Ao/eo — 1 in the heavy quark limit. The contents of these sum rules become more 
transparent in terms of the parton model. In parton language, the expectation values 
(O n ) of the operators composed of gluon fields can be expressed as Mellin transforms 
[12] of the gluon structure function of the light hadron, evaluated at the scale Q 2 = eg, 

(O n )= f dxx n - 2 g(x,Q 2 = e 2 ). (12) 

Defining now 

(13) 

we can reformulate eq. (11) to obtain 

/ dy y n - 2 <T* h (\ /y) = I(n) [ dx x n ~ 2 g(x,Q 2 = eg), (14) 
Jo Jo 



with I(n) given by 

16\ 2 r (n+ |) /16tt\ 1 



I(n) = 2n 3 / 2 



3 J r(n + 5) \3g 2 J m 



(15) 



Eq. (14) relates the $ — h cross section to the gluon structure function. To get a first 
idea of this relation, we neglect the n-dependence of I(n) compared to that of (O n ); 
then we conclude that 

cr$h(Xo/x) ~ g(x,Q 2 = Cq), (16) 

since all order Mellin transforms of these quantities are equal up to a constant. From eq. 
(16) it is clear that the energy dependence of the $— h cross section is entirely determined 
by the x— dependence of the gluon structure function. The small x behaviour of the 
structure function governs the high energy form of the cross section, and the hard tail of 
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the gluon structure function for x — > 1 determines the energy dependence of a<$>h close 
to the threshold. 

To obtain relation (16), we have neglected the n-dependence of the function I(n). 
Let us now try to find a more accurate solution of the sum rules (13). We are primarily 
interested in the energy region not very far from the inelastic threshold, i.e., 

(M h + e ) <A< 5GeV, (17) 

since we want to calculate in particular the absorption of $'s in confined hadronic 
matter. In such an environment, the constituents will be hadrons with momenta of at 
most a GeV or two. A usual hadron (tt, p, nucleon) of 5 GeV momentum, incident on 
a J/ip at rest, leads to y/s ~ 6 GeV, and this corresponds to A ~ 5 GeV. 

From what we learned above, the energy region corresponding to the range (17) 
will be determined by the gluon structure function at values of x not far from unity. 
There the x-dependence of g(x) can be well described by a power law 

g(x) = g 2 (fc + 1) (l-x) k , (18) 

where the function (18) is normalized so that the second moment (12) gives the fraction 
<72 of the light hadron momentum carried by gluons, < O2 >= #2 — 0.5. This suggests 
a solution of the type 

<T*h(y) = a(l-y) a , (19) 

where a and a are constants to be determined. Substituting (18) and (19) into the sum 
rule (13) and performing the integrations, we find 

r(q + i) = 1 2^ 2 g 2 \ / 16\ 2 r(n + f)r(fc + 2) 

T(n + a) \ m 2 Q )\3J \3g 2 J T(n + 5) T(k + n)' K ' 

We are interested in the region of low to moderate energies; this corresponds to relatively 
large x, to which higher moments are particularly sensitive. Hence for the range of n 
for which eq. (5) is valid [4], n < 8, the essential n-dependence is contained in the 
T-functions. For n > 4, eq. (20) can solved in closed form by using an appropriate 
approximation for the T-functions. We thus obtain 

I> + 1) n «-fc-5/2_ (21) 

T(A; + 2) v 1 

Hence to satisfy the sum rules (14), we need 

, 5 T(A; + 2) 
a = k + - a = const. —} (22) 

Therefore the solution of the sum rules (13) for moderate energies A takes the form* 
rr m torWa ( 16 Y f 1 **} 1 £(fc + 2) / A \ ^ 



* The same functional form for the energy dependence was also obtained in [2] 
through parton model arguments. 



To be specific, we now consider the J/V'-nucleon interaction. Setting k = 4 in accord 
with quark counting rules, using g 2 — 0.5 and expressing the strong coupling g 2 in terms 
of the binding energy eo (eq. 10), we then get from eq. (23) the energy dependence of 
the J/ipN total cross section 

/ \ \6.5 

<tj,4, n (\) ~2.5 mb x (l - -^J , (24) 

with A given by eq. (6) and Ao — (Mn + eo). This cross section rises very slowly from 
threshold, as shown in Fig. 2; for Pn ~ 5 GeV, it is around 0.1 mb, i.e., more than an 
order of magnitude below its asymptotic value. This clearly shows that hadronic matter 
at meaningful densities and for the volumes relevant for nuclear collisions is not efficient 
in breaking up quarkonium states. To illustrate this, we note that in a volume of the size 
of a uranium nucleus, the average path length for a J/ip is L ~ (3/4)(1.2 A 1 / 3 ) ~ 5.6 
fm. For matter of standard nuclear density uq = 0.17 fm -3 , a cross section (Jj/^n = 0.1 
mb leads to a survival probability S = exp(— ctj/^n L) ~ 0.99, and even at five 
times standard nuclear density, we still have a survival probability of more than 95%. - 
We should note that the high energy cross section of 2.5 mb in eq. (24) is calculated in 
the short-distance formalism of QCD and determined numerically by the values of m c 
and eo.* It agrees very well with the 2-3 mb obtained from photoproduction data via 
vector meson dominance [6]. 

For sufficiently heavy quarkonia, the calculation of the <E> — h cross section as pre- 
sented is firmly based on short-distance QCD. To check its applicability for charmo- 
nium, experimental tests are certainly important. This is indeed possible, since a com- 
pletely analogous derivation [4] provides the cross section of charm photoproduction, 
7 N — > cc X . Instead of eq. (24), we then get 

<r lN ^ c - c {v) ~ 1.2 fibx(l-^) 4 , (25) 

in terms of the scaling variable v = (s - M^)/2, with v = [(Mj/ip + M N ) 2 - M%\/2. 
The exponent in eq. (25) is again determined by the gluon structure function; it differs 
from that in eq. (24) because of the absence of the bound state Coulomb wave function 
in computing the coefficient function for cc photoproduction. In Fig. 3 we compare the 
result (25) to a summary of data [6-9]. Since eq. (25) contains the contributions of 
both open and hidden charm, we have reduced it in Fig. 3 by 25% to account for the 
absence of elastic shadow scattering in the data. It is seen to agree quite well with the 
experimental behaviour, which in particular shows clearly the small cross section in the 
threshold region. 

A similar test is in principle provided by J/i/'-photoproduction, 7./V — > J/ipN. The 
forward scattering amplitude of this process will give the total charm photoproduction 

* An alternative and in principle equivalent approach is to calculate the $ — N cross 
section by perturbative QCD [13]. The sum rule formalism used here has the advantage 
of providing the total cross section directly in terms of the largely model-independent 
scales rriQ, eo and the fraction g 2 . 
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cross section a 7 7v^ c5 , provided vector meson dominance holds and the real part of the 
amplitude can be neglected. Near threshold presumably neither is true, and so (Fig. 4) 
the result of our calculation falls there below the data [6], even though also here the 
cross section attains its high energy value very slowly. 

The calculation of the total $/i cross section as given above includes all possible 
quarkonium break-up mechanisms and becomes exact in the heavy quark limit. It is 
interesting, however, to consider also non-perturbative mechanisms made possible by 
the finiteness of the charmed quark mass. Here the "subthreshold" process [14] 

J/* + N -> A c + D (26) 

deserves particular attention, since it can occur even at very small energies. Apart 
from the perturbative contribution calculated above, the amplitude for reaction (26) 
can acquire a non-perturbative piece, corresponding to the tunneling of the c-quark 
from the charmonium potential well over a large distance Rt ~ A.Q^jp into the potential 
well of the nucleon. One can roughly estimate the probability of this process. The time 
r t required for the tunneling over the distance R t is r t ~ Aq^ d . On the other hand, 
the energy fluctuation needed to overcome the potential barrier of the height eo has 
according to the uncertainty relation a much shorter lifetime, of order Tf ~ e^" 1 . The 
probability for tunneling is thus expected to be 

W t ~ exp(-r t /r/) ~ exp(-e / Aqcd)- (27) 

This probability vanishes as expected in the heavy quark limit, since eo ~ tuq (see 
eq. (10)). For the J/ip, the probability (27) is already quite small (about 0.02), and 
hence the corresponding non-perturbative contribution is unlikely to affect the total 
quarkonium absorption cross section in any significant way. 

Since we have concluded that quarkonia cannot be broken up in hadronic matter, we 
should comment briefly on the implication of our results for hA collisions. In presently 
available data, the produced quarkonia (J/ip, T) have very high momenta in the 
rest frame of the target nucleus. Hence they are not yet physical resonances when 
they leave the nuclear medium; as an experimental consequence of this, J/ip and ip' , 
although very different in size, suffer the same nuclear suppression [15]. The observed 
distributions for these nascent quarkonia are in fact well described by a combination of 
quantum-mechanical coherence effects ("shadowing") [16] and colour octet interactions 
[17] of the QQ system in the nucleus. Thus there is so far no experimental information 
about the interaction of fully formed physical quarkonia in the confined medium of a 
heavy nucleus. To obtain such information, quarkonium production has to be measured 
at large negative xf values [17]. For a hadron beam incident on a nuclear target, this is 
difficult, since the slow decay dileptons are hard to measure. The relevant experiment 
is possible, however, for a nuclear beam incident on a hydrogen target, so that the 
study of quarkonia in nuclear matter becomes feasible with the advent of the Pb-beam 
at the CERN-SPS. Our analysis predicts no (i.e., less than 5%) suppression for the 
production of the basic (IS) J/ip in the range —1 < xf< — 0.4, in contrast to some 25% 
suppression obtained by using the asymptotic $/i cross section even in the threshold 
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region [18-20]. Since the mass of the ip' is just at the open charm threshold, it is much 
easier to dissociate and hence will suffer much stronger absorption. The same holds to 
a lesser extent also for the x states. 

In closing, we want to compare quarkonium interactions in confined and in decon- 
fined matter. In confined matter, the crucial gluon densities, g(x) ~ (1 — x) k , lead to 
(x) = l/(k + 2), which implies 

(P 9 ) = ^ (Ph) (28) 

for the average momentum of the gluons in terms of that of the hadron incident on the 
quarkonium. For hadronic matter of temperature T ~ 200 MeV, we thus get 

(p g ) ~ ~ 0.12 GeV (29) 

for the case of (massless) pions, with P n ~ 3T and k = 3, and 

(p g ) ~ \^2mT/it ~ 0.12 GeV (30) 

for nucleons, with Pn — 2^/2mI '/n and k = 4. It is clear that for gluons of such 
momenta it is difficult to break up quarkonium, i.e., to overcome the threshold to open 
charm. In previous absorption studies of charmonium suppression [21], it was suggested 
that this threshold is overcome in the collision of resonances (p, ...) with the J/ip, 
by making use of the difference between resonance and pion mass, M p — 2M^. Such 
arguments would be correct for the interaction between light and large hadrons; but in 
the interaction of the tiny quarkonia with light hadrons, the short-distance character 
due to the heavy quark mass and the large gap to open charm or beauty makes the 
gluon distribution within the light hadron the important feature, not its mass. 

In deconfined matter, on the other hand, we expect gluons in a medium of temper- 
ature T ~ 0.2 GeV to have an average momentum 

(p g ) ~ 3T ~ 0.6 GeV. (31) 

The impact of gluons with such a momentum spectrum can easily overcome the 0.7 GeV 
threshold to open charm and hence break up a J/ip into D's. The essential feature of 
deconfinement in this dynamical view of J/tj; suppression in a quark-gluon plasma [5] 
is thus the hardening of the gluon spectrum. 
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Figure Captions: 

Fig. 1: Deep inelastic scattering (a), heavy flavour photoproduction (b) and quarko- 
nium-nucleon interaction (c). 

Fig. 2: The energy dependence of the quarkonium-nucleon cross section, normalised to 
its asymptotic value; also shown is the dependence on the momentum Pjy of a 
nucleon incident on a J/tj; at rest. 

Fig. 3: The cross section of open charm photoproduction [6-9], compared to the pre- 
diction of the short distance QCD analysis. 

Fig. 4: The total cross section for charm photoproduction [6], compared to the predic- 
tion of the short distance QCD analysis, with the real part of the amplitude 
neglected. 
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